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We present a quantitative analysis of recent data on the kinetics of ATP hydrolysis, which has 
presented a puzzle regarding the load dependence of the Michaelis constant. Within the framework 
of coarse grained two-state ratchet models, our analysis not only explains the puzzling data, but 
provides a modified Michaelis law, which could be useful as a guide for future experiments. 
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Molecular motors or protein motors are the terms used 
to describe a highly specialized class of enzymes which 
transduce the energy excess in the chemical hydrolysis 
reaction of ATP (adenosinetriphosphate) into mechani- 
cal work. They are involved in several important cellular 
processes, ranging from the transport of material and 
vesicles, to cell mobility and cell division. Linear protein 
motors move along complex periodic and polar structures 
called filaments, obtained by the polymerization of a sin- 
gle monomer (actin) or dimer (microtubules) |^,^]. 

Several models have been proposed so far to explain the 
energy transduction process. In the earliest ones, known 
as cross-bridge models HB, the motor can exist in dif- 
ferent states (up to 5 or 6 [pi), within each of which the 
system reaches local thermodynamic equilibrium on time 
scales small compared to the exchange rates between the 
states. This hypothesis relies on the fact that the charac- 
teristic times of the motors (measured through transient 
response Q) are of the order of milliseconds, while ther- 
mal equilibrium on a length scale of about 10 nm occurs 
in 10-100 ns @J|. The latest models §-(j]| share some 
common features: i) need for asymmetry (polarity) in or- 
der to establish a certain preferred direction of motion; 
ii) chemical energy consumption as a source of mechan- 
ical work; iii) local thermodynamic equilibrium for each 
state. A particular promising approach is that due to 
Vale and Oosawa [ pL6|, using the ratchet concept intro- 
duced by Feynman [ |l7| . Since then, ratchet models have 
been intensively studied [p| JTc|jr^ ] . Non-equilibrium ther- 
modynamics as well as stochastic methods contribute to 
the theoretical background of this new field of research. 

The role of force in the reaction kinetics is still an open 
question, [l9|-|2l[|. In particular, recent experiments per- 
formed on kincsin (a processive motor) , [pll , seem to indi- 
cate that the kinetics of ATP hydrolysis can be described 
by the Michaelis-Menten mechanism. This simple mech- 
anism, introduced in 1913 to describe the process of en- 
zymatic catalysis J22], assumes that the catalytic reac- 



tion of a substrate S is divided into two processes. The 
enzyme E and the substrate first combine rapidly and 
reversibly to give an enzyme-substrate complex ES with 
no chemical change on the substrate. The chemical reac- 
tion occurs in a second step with a first-order rate con- 
stant k ca t, or turnover number. It is then simple to show 
that the reaction velocity r (or rate of S consumption) is 
given by the Michaelis law: 



[5'] 
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where r max is a saturation value and Km is called the 
Michaelis constant. For the ATP hydrolysis reaction, [S] 
is simply replaced by [ATP]. Furthermore it is reason- 
able that the velocity of a molecular motor should depend 
strictly on the rate of ATP consumption [Q], so that the 
velocity curve is given by a Michaelis law in terms of 
ATP concentration. This hypothesis has been tested in 
the experiment |2l) at least for high ATP concentrations. 
According to these experimental findings the Michaelis 
law can be written as: 
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where p is the periodicity of the filament, k cat is the 
turnover number, and e, the coupling ratio, gives an esti- 
mate of the average performance of one hydrolysis event, 
i.e. it is higher when the energy transduction process is 
more efficient. It is rather intuitive that the coupling 
ratio should depend on the applied force, and should 
decrease when the opposing force is increased, as the 
applied load strongly limits the maximum attained ve- 
locity on which the coupling ratio is linearly dependent. 
The rather surprising experimental finding pl| ] is that 
the Michaelis constant Km is no longer a constant, but 
increases with increasing applied load. This result is not 
intuitive and rather striking, since the Michaelis constant 
is just an equilibrium constant for the reaction leading to 
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the formation of the enzyme-substrate complex. Indeed 
reaction rates are assumed to be independent of the ex- 
ternal force as long as they do not contribute to a net mo- 
tion of the motor. Thus this experimental finding makes 
the picture more complicated than expected. It was sug- 
gested in [EH) that a possible explanation might be to 
insert the external force in the transition rates as done 
by Fisher and Kolomeisky pp| , but this cannot be done 
unambiguously, as stated by the authors themselves. 

Coarse-grained two-state ratchet models [Gl| can give 
further insight into this problem. In these models, the 
state of the motor is indicated by the index i, whereas x is 
the position of the protein center of mass along the track. 
The chemical reactions force the motor protein to switch 
from one state i at position x to another state j in the 
same position x, with a rate given by The motor 

moves under the influence of a potential Wi (x) chosen so 
as to reproduce the interaction between the filament and 
the motor head and with the same characteristics as the 
filament: polarity and periodicity. For this reason the 
potential is chosen to be asymmetric and periodic with 
period p for each state i ("ratchet" potential, see Fig. [I]). 

In a Fokkcr-Planck (FP) description, Pi(x, t) is the 
probability density for the particle to be in state i at 
position x at time t, while the probability current Ji(x, t) 
is defined as: 



Ji(x,t) 



A 

kT 



[Pi(x)d x Wi - P l (x)F ext + kTd x Pi(x)] 
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where A is a diffusion coefficient which is taken to be 
equal in both states, T is the temperature, k is the Boltz- 
mann constant and F ext is the external force. Since we 
are dealing with a simple two-state model, the chemical 
reaction cycle is compressed in a two step process: 



M + ATP M-ADP-P 



M + ADP + P ^ M ■ ADP ■ P 

02 (x) 



(4) 
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where M refers to the motor. This strongly resembles 
the Michaelis-Menten mechanism. The only difference, 
in this case, is that transition rates depend on the co- 
ordinate of the protein center of mass. We emphasize 
that all of the following discussion does not involve any 
hypothesis regarding enzymatic catalysis. 

The state M ■ ADP ■ P corresponds to one in which the 
motor head is detached from the fiber after binding and 
dissociating ATP. In this state the motor may move more 
or less freely upon the filament. This state will be called 
the "free" state, or state 2. All the other terms in the 
equations refer to a state in which the motor is attached 
to the filament and therefore its motion is strongly de- 
pendent on the motor-filament interaction. This is the 



so called "bound" state or state 1. In this case we are 
dealing essentially with a single-headed motor, although 
the mechanism for a two-head motor protein is similar. 
We assume detailed balance to hold for each chemical 
reaction: 
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where ai, a 2 (Pi, P 2 ) are the transition rates for the first 
(second) chemical reaction and A/i = ^atp — I^-adp — M-P 
is the difference in chemical potential or the chemical 
driving force. 

The FP equation describing the process may now be 
written as: 

d t Pi + d x J. t = J2 {"oiPj ~ UijPi) (8) 

where i = 1,2 for the two-state model and Wy(x) = 
oti(x) + Pi{x). 

The main advantage of this approach is the straight- 
forward manner in which external forces, if any, may be 
directly inserted in cq. (||) without any need for further 
assumptions (for a discussion of this point sec |20|). If 
the external force is independent of time, it is possible to 
look for a stationary solution. The velocity of the motor, 
ik defined as: 



(Ji(x) + J 2 (x)) dx, 



and the rate of ATP consumption: 



(ai(x)Pi(x) - a 2 (x)P 2 (x)) dx, 



(9) 
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where p is the periodicity of the filament, are identically 
zero when both A/i and F ext are zero ]18[| . A neces- 
sary condition for motion to occur |l8[ is that A/i ^ 0, 
i.e. detailed balance is violated, and that the potential is 
asymmetric. This implies that the chemical reaction of 
ATP hydrolysis is able to break detailed balance for the 
total transition rates. 

Using eq. (||) and (Q) only two of the four functions 
a\(x) , a 2 (x) , (3\(x) , and fi 2 (x), can be chosen arbitrarily 
once Wi and W 2 are fixed. Since the release of products 
is just a thermal process and does not involve chemical 
reactions (see eq. (||)), we assume it to be position inde- 
pendent as in [U, so that (3 2 (x) = u> = const. 

Transitions due to chemical reactions are usually cho- 
sen to be localized, i.e. they may take place only in corre- 
spondence to a certain motor position along the filament 
period and therefore they are not distributed over the 
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whole period. This corresponds to the "active site" con- 
cept in biology and in Jl8| the former hypothesis is shown 
to agree with experimental data. To take this effect into 
account we also define: 



a 2 (x) 



u> for p — 5 < x < p, 
otherwise. 
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with 6/p = 0.05. We use two kinds of models, differing 
only in the choice of the potential in state 2 [18|,||]. W\ is 
the standard ratchet potential shown in fig. IT In model 

(a) , we suppose state 2 to be strictly diffusive so that the 
potential W 2 {x) is flat (Fig.|l|). This corresponds to the 
picture of state 2 as a totally free state. In model (b), 
we suppose the filament to affect the protein movement, 
so that the potential in state 2 is essentially the same as 
W\ (x) , except for a uniform offset and a 5 times smaller 
amplitude of variation. We expect the motion in model 

(b) to be more constrained than in model (a). 

We analyzed the two models (a) and (b) seeking the 
non equilibrium stationary solutions P?(x) to the FP 
equations when applying different mechanical and chem- 
ical driving forces. This has been accomplished by nu- 
merically integrating the stationary FP eq. (^) starting 
from the equilibrium Boltzmann solutions and changing 
smoothly and alternatively the parameters Afi and F ext . 
For each stationary solution we calculated the velocity 
and the rate of ATP consumption. Fig. || shows the re- 
sults for the rate of ATP consumption. 

These results are plotted in terms of the ratio: 
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from mass-action law, using concentrations normalized 
with respect to their equilibrium value. Usually the 
Michaelis law is written in terms of ATP concentration, 
while in our results we derive it in terms of the parameter 
q defined above. We observe that experiments are usu- 
ally performed in conditions of high ATP concentrations, 
so that the ratio q can be safely thought to be constant 
during the time taken by the experiment. Since the only 
way of varying q is by adding ATP molecules to the so- 
lution, q or [ATP] may be interchanged in the Michaelis 
law. The rate of ATP consumption can be fitted by a 
Michaelis law in the form: 



q ■ T M AX 
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The value tmax is weakly dependent on the applied 
force. Our calculations show that the Michaelis-Menten 
law is followed with a rather impressive precision for 
1 ^ 10 3 , where 1/r is linear in l/q with regression coeffi- 
cients practically equal to 1. The fact that the Michaelis 
law is so strictly followed also when different forces are 
applied is a good indication that the simplified two-state 
model is still in accordance with well established results 
for enzymatic catalysis. 



Figure || shows the velocity as a function of q; at first 
sight these curves seem Michaelis-like, so that the correct 
equation for velocity should be eq. ([l3|) with tmax re- 
placed by vmax- In this case vmax strongly depends on 
the applied force, and it decreases when the applied force 
becomes more and more negative. Our results show that 
the velocity is linear with r, so that v = ar — (3 with a re- 
gression coefficient always > 0.9996 and thus the relation 
for velocity, using eq. (O), reads: 



otmax ■ 1 
K M + q 



0. 



(14) 



A pleasing feature is that, using this equation, there is 
only one Michaelis constant for both v and r. Another 
very interesting point is that the obtained positive val- 
ues of (3 indicate that ATP consumption is even possible 
under stall conditions (v = 0) |^3| . This "idling" rate of 
ATP consumption should be tested by motility assays. 

The Michaelis law is commonly thought to hold for 
very high ATP concentrations. The experimental re- 
sults pi| were obtained in this regime and fitted using 
eq. (|l4|) forcing [3 = 0. Table | shows the results ob- 
tained by fitting our numerical data using eq. (|lj) with 
(3 as a free parameter; these results are compared with 
the ones obtained in a fit with a standard Michaelis law, 
i.e. forcing /3 = 0. In the former case the fit is extremely 
good; indeed the sum of square residuals can be hundreds 
times less than in the latter case. While the maximum 
attained velocity is essentially independent of the fitting, 
the results for the Michaelis constant are quite different. 
The range of variability of Km with F ext is much wider 
using a standard Michaelis law; this is also confirmed by 
simple analytical arguments. Furthermore, in the case of 
model (b) for high loads it is even impossible to fit the 
data without introducing (3 in eq. (|l4|). Still, eq. (|^) does 
not allow for any inversion in the sign of velocity, nor for 
a stall condition, at variance with eq. (|l4|). 

In conclusion, our calculations on the two-state ratchet 
model clearly show that the rate of ATP consumption 
strictly follows a Michaelis law in the form (|13|); that 



the correct law for velocity is given by eq. (14), instead 
of eq. (0); that the observed fourfold increase of the 
Michaelis constant with applied load may be due to the 
use of eq. (||); that the model still accounts for an in- 
crease of the Michaelis constant smaller than 20% when 
forces vary by a factor 5 (in the range of pN), to be com- 
pared with the corresponding variation of about 250% if 
a standard Michaelis law, eq. (||), is used. 

We thank F. Cecconi and P. Van for useful discussions. 
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3 



[1] B. Alberts, D. Bray, J. Lewis, M.Raff, K. Roberts, and 
J.D. Watson. The Molecular Biology of the Cell. Garland, 
New York, 1994. 

[2] H.E. Huxley and W. Brown. J. Mol. Biol, 30:383-434, 
1967. 

[3] A.F. Huxley. Prog. Biophys., 7:255, 1957. 
[4] T.L. Hill. Prog. Biophys. Mol. Biol, 28:267, 1974. 
[5] N. Thomas and R.A. Thornhill. J. Phys. D, 31:253-266, 
1998. 

[6] A.F. Huxley and R.M. Simmons. Nature, 233:533-538, 
1971. 

[7] S. Leibler and D.A. Huse. J. Cell Biol, 121:1357-1368, 
1993. 

[8] A. Parmeggiani, F. Julicher, A. Ajdari, and J.Prost. 

Phys. Rev. E, 60:2127-2140, 1999. 
[9] M.O. Magnasco. Phys. Rev. Lett, 71:1477-1481, 1993. 
[10] J. Prost, J.F. Chauwin, L. Peliti, and A. Ajdari. Phys. 

Rev. Lett., 72:2652-2655, 1994. 
[11] C.R. Doering, W. Horsthemke, and J. Riordan. Phys. 

Rev. Lett, 72:2984-2987, 1994. 
[12] F. Julicher and J. Prost. Phys. Rev. Lett, 75:2618-2621, 
1995. 

[13] I. Derenyi and T. Vicsek. Phys. Rev. Lett, 75:374-377, 
1995. 

[14] R.D. Astumian. Science, 276:917-922, 1997. 
[15] A.B. Kolomeisky and B. Widom. J. Stat. Phys., 93:633- 
645, 1998. 

[16] R.D. Vale and F. Oosawa. Adv. Biophys., 26:97-134, 
1990. 

[17] R.P. Feynman, R.P. Leighton, and M. Sands. The Feyn- 
man Lectures on Physics. Addison-Wesley, Reading, 
MA, 1966. 

[18] F. Julicher, A. Ajdari, and J. Prost. Rev. Mod. Phys., 

69:1269-1281, 1997. 
[19] T. Duke and S. Leibler. Biophys. J., 71:1235-1247, 1996. 
[20] M. E. Fisher and A. B. Kolomeisky. Proc. Natl. Acad. 

Set. USA, 96:6597-6602, 1999. 
[21] K. Visscher, M. J. Schnitzer, and S.M. Block. Nature, 

400:184-189, 1999. 
[22] A. Fersht. Enzyme Structure and Mechanism. New York, 

Freeman, 1985. 
[23] We thank one of the referees for this observation. 



kT 




FIG. 1. Model (a): ratchet potential (Wi) and diffusive 
state 2. The parameters are chosen so that a/p — 0.1, 
U = WkT, U = 12kT,uj/D = 50. 
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FIG. 2. Results for the rate of ATP consumption for model 
(a) with lo/D = 50 and various forces (top to bottom curves: 
F ex t = 0.0, -1.0, -3.0, -5.0 in units of kT/p. 
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TABLE I. Michaelis constants, Km, and maximum veloc- 
ities, Vmax = a r max calculated from model (a) using eq. ( p^| ) 
(second and fourth columns) with /3 as a free parameter and 
the same equation with (3 = (third and fifth columns) cor- 
responding to the standard Michaelis law. In the latter case, 
the constants Km and v ma x are indicated as K' M and ul„„ 
respectively. Forces (first column) are measured in units of 
kT/p. 
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FIG. 3. Results for the sliding velocity for model (a) in 
units of p with the same values of the forces as in figj^. 
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